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ABSTRACT: 

We developed a biochemical kinetics approach to describe the repair of double strand 
breaks (DSB) produced by low LET radiation by modeling molecular events associated 
with the mechanisms of non-homologous end -joining (NHEJ). A system of coupled non- 
linear ordinary differential equations describes the induction of DSB and activation 
pathways for major NHEJ components including KU70/80, DNA-PKcs, and the Ligase IV- 
XRCC4 hetero-dimer. The autophosphorylation of DNA-PKcs and subsequent induction 
of yH2AX foci observed after ionizing radiation exposure were modeled. A two-step 
model of DNA-PKc S regulation of repair was developed with the initial step allowing 
access of other NHEJ components to breaks, and a second step limiting access to Ligase 
IV-XRCC4. Our model assumes that the transition from the first to second-step depends 
on DSB complexity, with a much slower-rate for complex DSB. The model faithfully 
reproduced several experimental data sets, including DSB rejoining as measured by 
pulsed-field electrophoresis (PFGE), quantification of the induction of yH2AX foci, and 
live cell imaging of the induction of Kuyo/so- Predictions are made for the behaviors of 
NHEJ components at low doses and dose-rates, where a steady-state is found at dose- 
rates of 0. 1 Gy/hr or lower. 
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INTRODUCTION 


A mechanistic description of the processing of DNA double strand breaks (DSB) is 
important for the understanding of ionizing radiation effects leading to cell death, 
mutation, genomic instability, and carcinogenesis. Mathematical models of DSB repair 
are important for the description of radiation modalities not accessible by experimental 
means and for their possible predictive capabilities. Past mathematical models of ionizing 
radiation induced DSB repair have largely relied on phenemological approaches, which 
did not consider specific molecular interactions involved in DSB repair (1-3). Previously, 
we had shown that a biochemical approach based on non-linear kinetics enjoys some 
special features in describing DSB repair, including the time delay caused by a DSB- 
repair enzyme intennediate complex (3). Application of biochemical kinetics models to 
describe molecular DSB repair experimental data is a goal of the present study. 

Non-homologous end joining (NHEJ) is the primary pathway for DSB repair in 
eukaryotic cells (4-6), and defects in NHEJ increase radiation sensitivity and the risk of 
carcinogenesis (7, 8). Many of the steps involved in NHEJ have been characterized 
experimentally, including the initial recognition of DSB’s by the KU 70/80 heterodimer, 
subsequent recruitment of the DNA dependent protein kinase catalytic subunit (DNA- 
PK CS ), and fonnation of the DNA dependent protein kinase (DNA-PK) (6-9). DNA-PKc S 
contains several serine-threonine residues which are auto-phosphorylated. Auto- 
phosphorylation of various subsets of these sites is thought to be important in regulating 
pathway choices between NHEJ or homologous recombination repair (HR) (6). In 
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addition to DNA-PKc S , a number of other proteins have been implicated as being 
important in either NHEJ or HR. For example, Artemis in conjunction with both ATM 
and DNA-PKcs has been suggested to function in DNA end-processing of specific, 
difficult to repair, IR-induced damages (10-13), and both Ligase IV and XRCC4 have 
been shown to be important in the ligation step of NHEJ (14, 15). 

DNA-PK CS is a member of the phosphoinositide-3-kinase-realted protein kinase (PIKK) 
family, which includes ataxia-telangiectasia mutated (ATM), and ataxia-telangiectasia 
and Rad3-related (ATR) proteins, and these proteins play a role in sensing DNA damage 
(5). In addition to damage sensing, both ATM and ATR have been shown to have distinct 
roles from DNA-PKcs, consisting of Gl/S, S and G2/M cell cycle checkpoints regulation 
(5, 16) and replication stress response (19), respectively. DNA-PKcs, ATM, and ATR 
share common features including a conserved carboxy-tenninal motif (5), and a reliance 
on upstream activators, Ku 7 o/so, the MRE-Rad50-Nbsl complex (MRN), and ATRIP, 
respectively. ATR is believed to function largely in S-phase, whereas DNA-PK CS and 
ATM have important roles throughout the cell cycle (16-19). The activation step of 
DNA-PK CS and ATM is rapid occurring from a few to about 30-minutes as observed in 
recent studies (5, 9, 18). Both activated proteins have been shown to lead to the 
phosphorylation of the histone variant H2AX in a chromatin region corresponding to 
about 2-Mbp around the DSB, with the phosphorylated form denoted yH2AX (20, 21). 
Total numbers of yH2AX foci have been shown to be fairly representative of the total 
number of DSB (20, 22, 23). In addition, a correlation between yH2AX foci loss and 
radiation sensitivity has been noted (23-25). 
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The plethora of experimental studies involving NHEJ repair should facilitate the 
development of mathematical models of these processes. In this paper we have developed 
a systems biology approach to NHEJ repair that can be used to make predictions for other 
radiation modalities, including extrapolations to low doses and dose-rates. Systems 
biology seeks to describe emergent properties of biological systems from the interactions 
of molecules acting in specific pathways (26). We use this approach to describe DSB 
rejoining curves, and the kinetics of formation and loss of yH2AX to gain insights into 
the kinetics of NHEJ repair pathway. A key component of a biochemical kinetics model 
is the role of DNA repair complex intermediates, which leads naturally to a non-linear 
kinetics description (3). We consider several intermediate complexes based on yH2AX 
radiation induced repair foci (RIRF) data, and DNA-PK CS experimental studies, and relate 
their description to pulsed-field gel electrophoresis (PFGE) DSB rejoining curves. 

Ionizing radiation produces DSBs that vary from simple to complex structures and are 
produced with equal proportions with increasing dose, and depend on radiation quality 
(27, 28). Clustered DNA damage sites are defined as two or more elemental lesions 
within one or two helical turns of DNA produces by a single radiation track (27, 29). 
Under this definition all DSB are clustered damages, however complex DBS’s are 
defined by the addition of other damage type such as base-damage, damaged ends, 
single-strand breaks near a DSB, or for two or more DSB in close proximity. Clustered 
non-DSB’s can lead to secondary DSB’s produced during damage processing (29, 30). 
For low LET radiation it has been estimated that 20-40% of the initial damage is complex 
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(27-29). Closely space multiple DSB could inhibit the attachment of repair proteins to 
other nearby DSB, and this possibility increases with the ionization power or linear 
energy transfer (LET) of radiation. We hypothesize that damage processing of complex 
DSB involve additional NHEJ factors including Artemis (10-13), MRN (19) and ATM 
proteins (10). 

METHODS 

DNA-PK Regulation and Repair Complexes 

We use the mass-action chemical kinetics approach to describe the binding of repair 
enzymes to DSB’s with several intermediate repair complexes leading to DNA rejoining: 
1) an initial complex bound by the KU 70/80 hetero-dimer, 2) Ku mediated DNA-PKc S 
binding, 3) The regulation of the DSB-DNA-PK cs complex through auto-phosphorylation 
by DNA-PK, and 4) a final repair complex involving the Ligase IV/XRCC4 heterodimer, 
denoted LilV. Figure 1 shows a schematic diagram of our model showing the sequence 
of proteins binding to the repair complex and the two activation steps considered; 
phosphorylation of DNA-PKc S and H2AX. The series of repair complexes are denoted Cj 
or with super-scripts P for auto-phosphorylation in complex, for e.g. Cf . The first 
complex (Ci) is formed by KU 70/80 binding to the DSB, and the second binding by DNA- 
PK CS to the first complex forming Cj etc., through to the final ligation step. Because these 
proteins are post-transcriptionally regulated the total number of enzymes in free-fonn or 
complex form is assumed to be conserved. 
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DSB are assumed to be induced per unit dose-rate with efficiency, a (Gy' 1 per cell). The 
KU 70/80 hetero-dimer is highly abundant and rapidly attaches to the DSB (denoted as Co) 
leading to the mass-action equation 


(i) = a ~r -*.[*“7.,«,][c 0 ] 

at at 


forming the repair complex, C/, which is followed quickly by DNA-PK C , binding 


(2) = ^[C 0 ][^ 70/80 ] -k 2 [DNAPKJ[C x ] 

at 


to form a second complex, C 2 . Equations (1) and (2) follow the convention that symbols 
within brackets define for a given molecular species, the time-dependent number of 
copies per cell. The repair complex, C 2 , is then modified by phosphorylation events that 
facilitate cleaning of the ends, signal transduction, and the translation of DNA-PKc S away 
from the ends of the break to allow ligation by LigaseIV/XRCC4 complex. Auto- 
phosphorylation of a cluster of residues on DNA-PK CS , denoted ABCDE, is expected to 
be a gate-keeper regulating access to the break by other repair proteins (6). The 
phosphorylation of a second cluster of residues on DNA-PKc S , denoted PQR, has been 
suggested to promote HR (6), whereas phosphorylation of the ABCDE cluster is thought 
to inhibit HR (6). The PQR cluster has been noted to partially facilitate dissociation of 
DNA-PK from the ends, however it is expected that other phosphorylation sites are 
needed for complete disassembly (6). To date, phosphorylation of K.U70/80 has not been 
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implicated as being critical for the actual repair of DSB (31) shown and thus will not be 
considered in the model. 

We consider a two-step model that depicts DNA-PK CS role in regulation of repair 
involving activation events controlled by auto-phosphorylation of DNA-PK CS . The exact 
nature of the auto-phosphorylation of DNA-PK CS is not known (6, 32, 33), 
autophosphorylation may occur in trans, where one DNA-PK CS molecule phosphorylates 
a second molecule on opposing sides of a DSB; a second-order reaction. Alternatively, it 
may occur in cis, by an intra-molecular mechanism; a first-order reaction. Both of these 
mechanisms may occur and could depend on the DSB end structure (32, 33). We have 
modeled the auto-phosphorylation of DNA-PK CS bound to the DSB ends as first-order for 
both steps in DNA-PK CS regulation of repair. We assume the second-step depends on the 
complexity of the DSB and may involve other proteins, including Artemis and other 
poorly defined repair proteins. Residual breaks are predicted at complex DSB sites 
through a competing first-order process that assumes not all complex DSB are 
successfully rejoined with the failure occurring before the transition to the ligation step of 
the reaction. The resulting equations are 

(3) d ^A = kl [DNAPKJ[C,\ - k Pi [C 2 \ 
dt 

d[C 2 P \ 
dt 


( 4 ) 


k Pi [C 2 ] - k Pi [C 2 P ]-k res [C 2 P ] 



(5) = k P \C 2 P ] - k,[LiIV][C 2 PP ] 


( 6 ) 


d[DSB r J 

dt 


= KJC/] 


The rates, k P 2 and /c; are assumed to depend on the complexity of the DSB, and k res set to 
zero for simple DSB. 

Finally, the last step in our model involves ligation of the ends by the Ligase IV/XRCC4 
complex, denoted UIV, and enzyme release given by 

(7) = k 3 [LiIV][C 2 PP ] - k Dc [C '] 

dt 

The Ligase IV/XRCC4 complex is also regulated by covalent modifications (13), but this 
observation is not currently treated in our model. 

yH2AX Foci Kinetics 

The histone variant H2AX is phosphorylated after DNA damage by each of the family of 
PIK3 phospho-proteins ATM, ATR, and DNA-PKc S (5). In the G1 phase of the cell cycle, 
ATM and DNA-PKc S phosphorylate H2AX with nearly equal efficiencies and in an 
overlapping manner (34). yH2AX foci appear at a distance from DSB corresponding to a 
region of 2 Mbp (20), and it is not known how many H2AX molecules are modified per 
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DSB or of the mechanism that leads to phosphorylation of a large number of H2AX 
molecules. We use Michaelis-Menten kinetics to describe the induction rate of H2AX by 
DNA-PKcs in its active forms as given by, 


( 8 ) 


d[yH2AX] 

dt 


k Pr [C Dm _ PK J[H2AX] 


+ \C DNA-PKcs . 


k Dy [yH2AX] 


where [Cdna-pkcs] is the sum of active forms of DNA-PKcs (C 2 P , C 2 PP , and C3). The 
mechanism of de-phosphorylation of yH2AX foci has not been well studied. We assume 
this step follows a simple first-order decay law in Eq. (8). 


Scaling Variables 

In order to simplify the model solutions, we introduce new scaled-variables by 
considering the conservation relations for the total concentration of a given protein, and 
noting the sum of all repair complexes is equal to the initial number of DSB. The new 
scaled variables are introduced using the definitions, 


H 


and 


[£,] + Etc,] = const.-, h t (t) 

j=' 



K, 


H,k , 


c ,(0 


[CJ 


which after substitution leads to the system of equations 
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( 9 ) 


dc 0 {t) 

dt 


-/fiC 0 (O(l-^i(O) 


a dD 
H l dt 


(10) ck j t l = -K oCx m-h 2 {t)) 

dt 

(11) = K 2 c x (t)(\-h 2 (t )) - k P c 2 {t ) 

dt 

dc P {t ) , . p, . 

(12) ~ = A: p c 2 (0 -{k Pi +k res )c 2 (t) 

( 13 ) d ° 2 ^ = k P c 2 (t ) - K 2 c 2 P (t)(\ - h 3 (t)) 

dt 

(14) = k c ' (t) 

dt 

(15) ck j t l = ^(0(1 -^(0) -k D c 3 (t) 

dt 

All rate-constants are assumed to be independent of the type of initial DSB, except for kp 2 
and K 3 , which are given distinct values for simple and complex DSB, respectively. For 
the solutions in terms of the scaled-variables, only the value of Hi enters as all other II, 
are combined with the A:, to form the rate parameters, k, which are in units of Hr' 1 . The 
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functions h,{ t) include contributions from repair complexes involving both simple and 
complex DSB’s. 

The histone variant, H2AX content varies with cell lineage representing from 2 to 10% of 
all nucleosomes and there are about 2.0 xlO 6 H2AX molecules per cell (20). We reasoned 
that it was more useful to model the kinetics of the number of yH2AX foci fonned, rather 
than the number of activated molecules. For foci counting experiments the number of foci 
is limited by the model dependent initial number of DSB per cell. For low LET radiation 
the probability of more than one DSB within the spatial region of foci is small, however 
for high LET other considerations will be needed to be taken into account (Cucinotta et 
al., in preparation). Assuming [yH2AX]+[H2AX] = constant, and denoting y(t) as the 
time dependent number of foci leads to 


(16) 


dyit) 

dt 


^P y C DNA-PKcs (0(1 -K0) 


K,, + C , 


M) 


-k D y(t) 


where k p = k Py / // . For comparison to DSB rejoining kinetics in an acute irradiation 
measured using PFGE, the number of DSB’s remaining is given by 

(17) DSB remaining (t ) = I/ l [Y J c l (t) + c r Jl)] 

j = o 

For comparison to experimental data on relative KU 70/80 induction, which includes KU 70/80 
in various DSB repair complexes the following sum is used 
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(18) Q, 0 

;'= i 

The system of equations formulated above to represent NHEJ are non-linear ordinary 
differential equations, described as stiff equations describing equations were the values 
for the various parameters, kj or k, vary over several orders of magnitude. These 
equations were solved numerically using the method of backward difference 
approximates. We note that the factors ‘ /-hf in our scaled equation have values close to 
unity at low doses where the initial number of DSB is «Hi. 


RESULTS AND DISCUSSION 


Our kinetics model of NHEJ consists of a system of 8 coupled non-linear ordinary 
differential equations for each class of DSB (simple and complex). This system of 
equations describe major components in the NHEJ repair pathway and the 
phosphorylation of H2AX by DNA-PK CS . Values for rate-constant were determined by 
comparing to experimental data with cell lineage specific values estimated for rate- 
constants and other parameters are listed in Table 1 . Our scaling approach results in a 
significant reduction in parameter space since it avoids the need to estimate values for the 
total cellular concentration of K.U 70 / 80 , DNA-PK CS , LilV, and XRCC4, which are 
effectively replaced by a single constant, Hi The value of Hi can be interpreted as the 
total number of copies of Ku 7 o/8o- However in the model other constants, Hj, could be 
used as the scaling variable, and we prefer to interpret the value of Hi as the total 


13 



number of DNA-repair complexes that could occur in a cell (3). We have fixed this value 
at a large number (7//=3000), to ensure that the shape of the DSB rejoining curve is 
largely independent of dose, over the range from 1 to 40 Gy. To reduce the number of 
variable parameters, we fixed the peak time of the [CJ complex, corresponding to the 
binding of Kuyo/soto DSB, at about 1 min post-irradiation (35), and of the [CJ complex, 
corresponding to the binding of DNA-PK^ complex, at about 3 min for all cell lineages 
considered using the values for Kj and k 2 as listed in Table 1. In-vitro assays provide 
insights into rates for DNA-PKcs activation occurring over times up to 30 min under 
different conditions (32, 36, 37). The remaining parameters are determined in a cell- 
lineage specific manner by comparing the model solutions to data for DSB rejoining, and 
the induction and loss of yH2AX foci. 

Figures 2 illustrates the model predictions for the time evolution for the sequence of 
repair complexes formed at an acute dose of 1 Gy for simple and complex damage 
processing. We compared the model prediction to recent results using live cell imaging of 
DSB induced by a near infrared laser (NIR) of KU 70/80 (35) in Figure 3. Our result using 
eq. (18), which represents the sum over different repair complexes containing the KU70/80, 
is in excellent agreement with the live cell imaging observations where DSB’s are 
induced by a near infrared laser (NIR) (35). There will be differences in the initial 
number and types of breaks, perhaps due to higher DSB induction by NIR as compared to 
that of X-rays at 1 Gy, however the agreement found over the first few hours of repair 
lends support to the values for the rate-constants chosen. 
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We also compared our model to the rejoining kinetics detennined by PFGE, which is 
available in the literature for X-rays. DSB rejoining kinetics measured by PFGE are made 
at high dose (>10 Gy) and must be corrected for the presence of heat-labile sites (29, 38, 
39), which account for up to 50% of the fragment yields at early post- irradiation times 
(within 30 mins post-irradiation). To avoid the contribution of heat-labile sites, we 
compared PFGE data analyzed using the cold lysis method developed by Rydberg (38) to 
the data of Sternerlow et al. (39) for the GM5758 diploid fibroblast cells as shown in 
Figures 4. We have used experimentally determined values from Gulston et al. (29) for 
HF19 cells for the value of a, the total number of DSB per Gy, of 25 and 16, 
respectively as initial conditions, and assume that 20% of the initial breaks repair with 
additional processing steps between the transition from Ci P to C 2 PP and hence slower 
kinetic parameters. At moderate doses (<5 Gy) the model predicts a lack of rejoining in 
the first few minutes post-irradiation as the multiple steps in NHEJ proceed. However, 
there may be some DSB rejoined by direct ligation independent of DNA-PK CS (35) 
leading to a faster component at early times, beyond those contributed by heat-labile 
sites. 

Track structure calculations provide some estimates of the fraction of simple and 
complex DSB lesions. However, the bevy of mechanisms that would be available to 
repair the differential spectrum of DSB produced by ionizing radiation are not well 
understood and may utilize additional factors amongst these being Artemis (10-13), ATM 
(10), MRN (19, 40), Werner syndrome (41) proteins, and perhaps components in the 
nucleotide or base- excision repair pathways. In our model we assume just two average 
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components corresponding to so-called broad categories namely simple and complex 


DSB. 


The fraction of residual breaks is easily modeled when the complex DSB are considered, 
if one assumes a first-order process results near the end of the cascade described above. 
We used a first-order rate-constant for residual break formation of 0.05 If 1 assuming a 
small fraction corresponding of the initial complex DSB are remain un-repaired at the C 2 
complex and lead to residual DSB. The model presented here thus provides a framework 
to describe the dependence of residual breaks on radiation quality, dose-rate, and post- 
irradiation time. 

We compared our model to data for the time courses and dose-response for yH2AX foci. 
Leatherbarrow et al. (23) using confocal microscopy precisely measured the number of 
yH2AX foci in V79 and HF19 cells. We find good agreement with their results as shown 
in Figure 5. Comparisons of the number of yH2AX foci at 0.5 and 4 hr post-irradiation 
made by Short et al. (42) are shown in Figure 6. The model calculation shows a linear 
response at 0.5 hr post-irradiation. There is a concomitant induction of yH2AX from 
active ATM (34) monomers, which has not been studied in the current model. ATM and 
DNA-PK CS are expected to induce these foci with nearly equal efficiency (34). 

The understanding of dose-rate effects is an important consideration in radiation 
protection (43, 44). Since the processing of DSB after radiation is a detenninant in 
mutation, chromosome aberrations, and carcinogenesis, we studied the induction of 
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various NHEJ components as a function of variable dose-rates and doses. Steady-state 
solutions for the systems equations can be found in closed form and compared with 
numerical solutions, and dose-rates where the steady-state, with foci counts independent 
of dose-rate, are obtained identified. The results of Figure 7 predict that the number of 
DSB repair complexes per cell becomes independent of dose-rate below about 0.1 Gy/hr. 
These observations can be tested with experiments. Also, for model building our 
description of NHEJ can be used as a starting point for mechanistic models of mutation 
and chromosomal aberrations. In addition, studying dose-rate dependencies for repair 
should be informative in understanding dose-rate effects for these other endpoints. 

In summary, we have synthesized a large number of experimental observations into a 
biochemical kinetics model of the NHEJ repair pathway. The model is based on the 
current mechanistic understanding of molecular binding and kinase activity of major 
NHEJ components that have been described experimentally and is capable of describing 
the time-courses, and dose and dose-rate dependencies for major NHEJ components, the 
induction of yH2AX foci upon activation of DNA-PK CS through auto-phosphorylation, 
and DSB rejoining curves as measured by PFGE. The model presented here can be 
modified as understanding on molecular mechanisms of NHEJ repair is obtained. The 
ability to describe the kinetics of DSB induction and repair and the various associated 
protein complexes will support models of chromosomal aberrations as a function of 
radiation quality, when descriptions of DSB complexity and spatial dependence of initial 
DSB are coupled to the present model (27, 28, 45). We plan on extending our work to 
include theoretical descriptions of the fractions of simple and complex initial DSB for 
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high LET radiation, and the resulting changes in DSB repair and foci kinetics and to 
include the description of the ATM signaling pathway in our model. 
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Table 1. Values of Rate-constants and other parameters in the Biochemical Model 

a),b) 


Rate Constant 

V79 Cells 

HF19 cells 

T98G 

a, Gy' 1 

16 

25 

25 

K 3 .hr " 1 

8 (0.5) 

8 (0.5) 

8 (0.5) 

kpi.hr ' 1 

10 

10 

10 

kp 2 .hr ' 1 

10(0.5) 

10(0.5) 

10(0.5) 

Kp y . copy ' 1 hr ' 1 

1000 

900 

1000 

k Dc , hr " 1 

4 

4 

2 

kp)y. hr ' 1 

2 

2 

0.75 

Km 

0.5 

0.5 

0.5 

kres? hr 

0 (0.05) 

0 (0.05) 

0 (0.05) 


a) Values of Ki , K 2 , and Hi are set at 100 hr' 1 , 100 hr 1 , and 3000 per cell for all cell 
lineages considered. Values chosen correspond to a peak for a DSB-KU70/80 
complex at about 1 min post-irradiation ( 35 ), and assuming peak of DSB-KU70/80- 
DNA-PK CS complex at about 3 min post-irradiation. 

b) The initial number of breaks per Gy (a) determined from experiments of Gulston 
et al. (34). We use the same values for GM5738 cells as HF19. 
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Figure Captions: 

Fig. 1. Schematic of biochemical kinetics model of DSB repair by NHEJ with induction 
of yH2AX by DNA-PK CS . The key components of the model and associated rate-constants 
are show. Not illustrated is the degradation of the [C 3 ] complex after the ligation step nor 
distinction between simple and complex DSB including the fonnation of residual DSB 
when complex initial DSB do not proceed to the [C 2 PP ] complex. 

Fig. 2. Model calculations of time course for sequence of DNA repair complexes in 
NHEJ pathway, and DSB rejoining curve (non-complex only) for 1-Gy gamma-ray 
exposures in normal human diploid fibroblast cells. 

Fig. 3. Model calculations for the time-course of Ku70/80 hetero-dimers in complex with 
DSB and other NHEJ components compared to live cell imaging data from ref. (35) for 
EGFP-KU 80 induction after irradiation of CHO cells with a near infrared laser. 

Fig. 4. Comparisons of model calculations to DSB rejoining determined by PFGE method 
for GM5758 (human diploid fibroblast cells) at 40 Gy (37). The solid line shows the 
contributions from simple and complex DSB. For comparison we show calculations of 
the rejoining curves for simple and complex DSB in our model normalize to unity as dash 
and dotted lines, respectively. 

Fig. 5. Comparisons of model calculations to measurements of yH2AX foci by 
Leatherbarrow et al. (23) with solid line for total (simple and complex DSB) induced 
yH2AX foci, dash line foci from simple DSB alone, and dotted line the number of DSB 
remaining. Symbols with error bars are the experimental results (23). 

a) HF 19 Cells at 1 Gy 

b) V79 Cells at 1 Gy 

Fig. 6. Comparisons of model calculations for dose-response for yH2AX foci at 0.5 
(closed circles) and 4 hrs (open triangles) post-irradiation data of Short et al. (34). Model 
calculations shown are solid line at 0.5 hr, dash line at 4 hr post- irradiation. 

Fig. 7. Predictions for the number of DSB repair complexes versus dose for various dose- 
rates in human fibroblast cells. 
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